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Abstract 

We construct an equation-of-state for nuclear matter based on the chiral Lagrangian. 
The relevant scales are discussed and an effective chiral power expansion scheme, 
which is constructed to work around the nuclear saturation density, is presented. 
A realistic equation-of-state is obtained by adjusting one free parameter, when the 
leading and subleading terms in the expansion are included. The saturation mecha- 
nism is due to correlations induced by the one-pion-exchange interaction. Further- 
more, we find a substantial deviation from the Fermi-gas estimate of the quark 
condensate in nuclear matter already at the saturation density. 



For more than four decades, one of the basic problems of nuclear physics has 
been to find a microscopic understanding of the nuclear equation-of-state in 
terms of the free nucieon-nucleon interaction [1-3]. The present status is that a 
quantitative description of nuclear matter can be achieved in a non-relativistic 
approach only by invoking a three-body force [5], or alternatively with the 
lowest-order term in the relativistic Brueckner approach [4]. A drawback of 
the latter, is the lack of a systematic expansion scheme. Thus, it has so far not 
been possible to obtain a reliable estimate for the corrections to the leading 
term. Such calculations are based on the traditional semi-phenomenological 
approach to the nucieon-nucleon interaction, represented by e.g. the Urbana- 
Argonne [6] and Bonn [7] potentials. 

During the last few years, a novel approach to the nucieon-nucleon interac- 
tion, based on chiral perturbation theory (xPT), has started to emerge. A 
key element in is the power counting, which relies on a separation of 
scales. This allows one to organize a systematic approximation scheme. In the 
nuclear many-body problem, new scales enter, which requires a modification 
of the expansion scheme developed for the nucieon-nucleon interaction in free 
space. In this letter we identify the new scales and present an effective chi- 
ral perturbation theory, appropriate for the nuclear many-body problem. We 
employ this scheme to construct a nuclear equation-of-state and compute the 
quark condensate in nuclear matter. 
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The key element of any microscopic theory for the nuclear equation-of-state 
is the elementary nucleon-nucleon interaction. In the context of chiral per- 
turbation theory this problem was first addressed by Weinberg who proposed 
to derive a chiral nucleon-nucleon potential in time ordered perturbation the- 
ory [8,9]. An alternative scheme, with chiral power counting rules applied 
directly to the nucleon-nucleon scattering amplitude, was subsequently pro- 
posed [10-12]. This approach relies on the crucial observation that the chiral 
power counting rules can be generalized for 2-nucleon reducible diagrams. 
Non-perturbative effects like the pseudo bound state pole in the 1 S , channel 
are generated by properly renormalized local two-nucleon vertices, that carry 
an anomalous chiral power Q^ 1 . 

The counting rules for the vacuum nucleon-nucleon scattering amplitude pro- 
vide a suitable starting point for the construction of chiral power counting rules 
for the nuclear matter problem [13]. The density expansion, or equivalently the 
multiple scattering expansion, is readily combined with the chiral expansion, 
once one identifies the Fermi momentum /jp ~ Q as a further small scale. This 
identification leads to a systematic partial resummation of the density expan- 
sion, where one avoids the expansion in large ratios like kp/m-,, but exploits 
the chiral mass gap and expands in small ratios like &f/to x and m n /m x with 
m x ~ 4 7r f n ~ 1 GeV. 

In the generalized chiral power expansion scheme [10-12], the nucleon-nucleon 
scattering amplitude is expanded in the small ratio Q/A L , but not in the large 
one Q/As, where A L ~ m x and A s is a small scale, like ^m,N £d where — 2 
MeV is the binding energy of the deuteron. This leads to an expansion of the 
scattering amplitude of the generic form 



where the small scale Q is identified with the nucleon momentum or the mass 
of the pion m n . An immediate consequence of the chiral expansion schemes, is 
that the long-range part of the one-pion-exchange interaction can be treated 
perturbatively. The short-ranged or large-momentum part is, as described be- 
low, effectively included in the local part of the chiral Lagrangian, which is 
treated non-perturbatively. 

When this scheme is applied to the nuclear many-body problem one finds that 
the pion dynamics remains perturbative in the sense discussed above, but the 
local two-nucleon interaction requires extensive resummations. This is an im- 
mediate consequence of the chiral power given to the renormalized two-nucleon 
coupling g R ~ Q^ 1 . Thus, in nuclear matter at small density, say p ~ 0.01 
fm~ 3 , it is not sufficient to sum only the particle-particle ladder diagrams of 
the local two-nucleon interaction, as is done in lowest-order Brueckner calcu- 
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lations. We find that one must sum also the particle-hole and hole-hole ladder 
diagrams including all interference terms with self-consistently dressed nu- 
cleon propagators. According to the generalized counting rules the pions can 
then be evaluated perturbatively on top of the parquet resummation for the 
local two-nucleon interaction, described aboveQ. In terms of vacuum scales 
this leads to an expansion of the energy per particle, E(kp), of the form 



E{k F )=Y J E, 

n 

The expansion coefficients E n are complicated, presently unknown, functions 
of the Fermi momentum kp. However, by employing appropriate resumma- 
tions, these functions can be computed from the free-space chiral Lagrangian [14]. 

Since the typical small scale As — 50 MeV is much smaller than the Fermi 
momentum at the saturation density, kp ~ 265 MeV, one may expand the 
coefficients E n around h$ in the following manner 
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Note that we do not expand in the ratio m^/kp. Clearly this scheme is con- 
structed to work around nuclear saturation density but will fail at very small 
densities. However, also other approaches, like the Walecka model [17] and the 
Brueckner scheme, are not reliable at very low densities. 

Physically the expansion (3) may be interpreted in the following way. The 
parquet resummation of the local nucleon-nucleon interaction gives rise to a 
density dependence governed by the scale As, while the relevant scale for the 
pion dynamics is m n . If As/kp^ is so small, that the parquet resummation is 
close to its high-density limit at normal nuclear matter density, the expansion 
(3) converges rapidly, and we can, as a first approximation, retain only the 
leading term in this expansion. In this case the density dependence of the 
coefficients of the expansion (2) is, to a very good approximation, given by 
the pion dynamics. 

A systematic derivation of the expansion (2) and (3), applying suitable resum- 
mation techniques, will be presented elsewhere [14]. In this work we pursue a 



1 The expectation that the nuclear equation-of-state can be computed microscopi- 
cally by summing the parquet diagrams for the local interaction and then including 
pions perturbatively may be too optimistic. A refined formulation, which incorpo- 
rates the relevant subthreshold singularities of the vacuum nucleon-nucleon scatter- 
ing amplitude (see [11]), may require a more involved treatment of pionic effects. 
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less microscopic approach. We assume that the coefficients E n (kF/m n , hp /As) 
can be computed in an effective theory whose unknown parameters are ad- 
justed to the saturation properties of nuclear matter. We employ an effective 
Lagrangian density, where the nucleons interact through s-wave contact in- 
teractions and through the exchange of pions. Here we do not include local 
p-wave nor three-body interaction terms since they do not contribute to lead- 
ing and subleading order in the chiral expansion. The interaction part of the 
Lagrangian is given by 

bnt(k F ) = -gj-tf (<7-V) (tT-tW 

^ J IT 

+ ^ (m(k r ) + \A) (n'n) i (ivt JV ) 2 (4) 

where N is the two component nucleon spinor field and P x f T is the projection 
operator for a two-nucleon state with spin S and isospin T. Note that in the 
contact interaction the spatial coordinates of the two particles are identical. 
The terms proportional to g\/ 4 are counter terms, which cancel the local, high 
momentum, piece of the one-pion-exchange interaction. Below we discuss the 
role of the counter terms in more detail. 

Rather than evaluating the coefficients of the chiral density expansion (2,3), 
we compute the energy directly. The coefficients E n can be extracted from the 
resulting expressions. In Fig. 1 we show all diagrams of chiral order Q 3 and 
Q A . They correspond to the leading and subleading interaction contributions 
to (2). The dashed line is the pion propagator. We split the non-interacting 
nucleon propagator into a free-space part 

S N (po, P ) = ^W, ( 5 ) 

p -p 2 /(2m N ) +te 

denoted by a directed solid line, and a density-dependent part 

AS N (p ,p) = 2iri5( Po -p 2 /(2m N ))e(k F -p 2 ) (6) 

represented by a line with a 'cross'. Thus, depending on the diagram, such a 
line corresponds to a hole line or the Pauli blocking of a particle line. This 
separation is useful for the expansion scheme we adopt. 

The pion dynamics which is computed explicitly in our scheme. If properly 
renormalized, it is perturbative like in the vacuum case. The non-perturbative 
short-range physics is subsumed in the local interaction. In our scheme, the 
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power counting is much simpler than in a fully microscopic one, where one 
would have to sum diagrams involving the small scales As and kp to all orders. 
This infinite set of diagrams is already included in our interaction constants 
go(kp) and g\{kp). Consequently, in order to avoid double counting, diagrams 
with two or more adjacent interaction vertices go or g\ are forbidden. Therefore 
we introduce two types of local 2-nucleon vertices in Fig. 1. The filled circle 
represents the full zero-range vertex of (4) proportional to go,i + #1/4 while 
the open circle corresponds to only the counter term proportional to g\/^- 
Clearly there are no diagrams of the type e-f in Fig. 1 with the full contact 
interaction (filled circles) at both vertices. 

Note, that in our counting g and g\ are of order Q° since the non-perturbative 
structures like the deuteron, which give rise to the anomalously large ampli- 
tude ~ Q^ 1 in vacuum, are dissolved at densities far below the saturation 
density. We neglect the density dependence of the effective coupling constants 
go and g±. This corresponds to retaining only the leading term in (3). We can 
test whether subleading terms are required, by allowing the coupling constants 
to be density dependent. 

The contribution of the first diagram in Fig. 1 to the ground state energy 
density e is proportional to (go(kp) + k F . Since the effective vertices 

g (kp), gi(kp) ~ Q° carry chiral power zero, this diagram is of chiral order 
Q 6 . The contribution to the energy per particle E{k F ) = e(k F )/p is of order 
Q 3 , since p = 2k F /(3ir 2 ). In this paper, the chiral power of a diagram is 
defined by its contribution to the energy per particle. The one pion exchange 
contribution (Fig. 1 b), is also of chiral order Q 3 since it is proportional to k 3 F 
times a dimension less function of kp/m^. 

The leading interaction contributions to the energy per particle (Fig. 1 a,b) can 
be expressed in terms of a momentum dependent effective scattering amplitude 
Teffip): 



E T {k F ) = -^- 2 J x 2 dx{\-xf{2 + x)T eff {kpx) (7) 



where 



In (8) the cancellation between the counter terms in the local effective inter- 
action (4) and the high-momentum part from the one-pion-exchange interac- 
tion is evident. As an example for the link between the Lagrangian (4) and 
the chiral density expansion, we consider the coefficient E 3 [k F /m n ,k F ^ / A s ] 
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in (3). The contribution of the contact interaction is given by E^/K\ = 
— 2(g + g 1 )/(47r/, r ) 2 . We note that the small scale Ag is hidden in the coupling 
constants g and gi. 

Next we consider the diagrams of sub leading order, shown in Fig. 1. These 
diagrams fall into two distinct classes: those with two and three crosses. We do 
not show the diagrams with four crosses since their real parts vanish and their 
imaginary parts cancel that of the other diagrams so that, to a given order, the 
resulting energy is real. Furthermore, we do not include vacuum-polarization 
diagrams and vacuum vertex renormalization diagrams. The former vanish 
in a non-relativistic scheme^, while the latter are suppressed by two chiral 
powers Q 2 . 

We note that after a proper renormalization, diagrams with only one cross are 
implicitly included in the kinetic energy term 

£kin(M = ^— ■ (9) 

10 m N 



The contribution of diagrams d), f), h), j) and 1) with two crosses can again 
be expressed in terms of an effective scattering amplitude T e ff(p) with 



T^ f \p) = Q 9 -i(go + 9i + g -i\ n,il(,» 



/A 4 
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effW 2/ 4 4?r I l 6pm 3( m 2 +2p 2) I ^ f m 2 



pi p 2 \ \ ( p 



+ arctan — 3 + 4 — - — arctan 

*P ^~~(^ o • P \\ 3 9A 



2m 2 V mj) 2/4 



log 1 - 2 i ±- - ^ m% (l(p) - 2 1(i m n )) . (10) 



2 The non-relativistic scheme presented here is identical with the non-relativistic 
reduction of a relativistic scheme order by order in the chiral expansion [11]. 
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a) 




Fig. 1. Leading and subleading contributions to the energy of nuclear matter. The 
Feyman rules are explained in the text. 

The loop integral I(p) is divergent and must be regularized: 



I{p) = 



J 



l7r 2 rriM J I 2 — p 2 
o 



i e lQnm N \n \A 2 )J 



Here A is a cutoff parameter. 

In (10) the appropriate symmetry factors are included. We note that the dia- 
grams d), f), h), j) and 1) in Fig. 1 are all divergent. However, the sum of all 
diagrams is finite once we include a further counter term in the spin triplet 
channel: 



9o^ 9o~ 



m N A g\ 
4vr 2 f 2 



(12) 



This counter term, which is formally of order Q, cancels the divergence due 
to the second-order tensor interaction. Collecting all real terms from (8) and 
(10) we find 



Teff(p) = ^(go + 9i) + l^ 
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+ ^i{9o + 9i) 1R , 2 — arctan — 

m n m N ( ml _ 3 &rctan /2p 
/ 2 167r/2 \^ m 2 + 4p2 2 p Vm^ 



+ 



3g\ m n m N ( p_ /2p^ _ 1 

/ 2 167T/2 \ m w 



+ m * +8 ™y + y f arctan (p_ ( 3+4 p 2 



8pm ff (m 2 + 2p 2 ) \ 



m, 



2 



- arctm (^))j- (13) 

The corresponding contribution to the energy per particle is again given by 
Eq. (7). 

We now turn to the remaining diagrams c), e), g), i) and k). Their contribution 
to the energy per particle reads 



m N k% (3 2 / ^\ 
- | ^ (<to + gi + f ) Ji(*f) + \g\ Mk F ) - \g\ Hk F )\ (14) 



where 
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and 



A(q, x q ) = ^/fc F -q 2 (l-a;2) - |q| x q 
l 2 

W) 



2 

l 2 m 2 



F 1 (l) = l -- , -!±]og(l* + m£) 
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m = \ - 2 { ™; ml) - ml log (l 2 + ml) . (16) 

The integral Jo can be performed analytically with the result Jo = 16 (11 — 
ln4)/35. 

The resulting equation-of-state is given by E(kp) = Ej, in + Ep + Ep. We note 
that in symmetric nuclear matter there is effectively only one free parameter, 
namely the effective coupling constant g = go + gi, since we neglect the density 
dependence of go and g±, as discussed above. In Fig. 2 we show the energy per 
nucleon for g = 2.8, 3.0, 3.2, 3.4 and 3.6. We emphasize that the coupling 
functions (?o(&f), A , i(^f) are to be determined from the nuclear equation-of- 
state. Thus, by allowing <7o(&f) and gi(kp) to be arbitrary functions of fop 
we could have trivially obtained a realistic equation-of-state. However there 
is a strong consistency constraint: according to our scale argument discussed 
above, the density dependence of the coupling functions must be weak for 
kp larger than the small scales, which formally have been integrated out. 
Thus, our scheme would have to be rejected, had we found a strong density 
dependence of go and g±. In fact, the density independent set of parameters 
g ~ 3.23, gA — 1-26, m % ~ 140 MeV and f n ~ 93 MeV yields an excellent 
equation-of-state. We obtain saturation at the Fermi momentum kp^ ~ 265 
MeV, which corresponds to a density po ~ 0.16 fm -3 , in agreement with 
the empirical value. Furthermore, the empirical binding energy of 16 MeV is 
reproduced and we obtain an incompressibility of k ~ 218 MeV, compatible 
with the empirical value (210 ± 30) MeV of ref. [16]. 

It is interesting to explore the convergence properties of the expansion (2). At 
the saturation density, the leading term of chiral order Q 2 (the kinetic energy) 
contributes 22.5 MeV to the energy per particle E(kp), the terms of order 
Q 3 (Fig. 1 a-b) contribute -93.5 MeV and the terms of order Q 4 (Fig. 1 c-1) 
+55.0 MeV. We observe a partial cancellation between the Q 3 and Q 4 terms. 
However, based on these terms alone, one can draw only qualitative conclusions 
on the convergence properties of the expansion. Quantitative results have to 
await the evaluation of terms of higher order in Q. As one finds e.g. for the 
hole-line expansion in Brueckner theory [2], it may well be that the convergence 
is in fact much better than suggested by the first few (fairly large) terms. 

The saturation mechanism of our model is quite different from that of popular 
effective models, like the Skyrme and Walecka models. In the Skyrme model 
the saturation is due to the density and velocity dependence of the zero-range 
effective interaction [16], while in the Walecka model it is a relativistic ef- 
fect [17]. In neither of these models are pionic degrees of freedom explicitly 
included. In our model, on the other hand, saturation is obviously due to pion 
dynamics, since for g^ — only the diagram in Fig. la survives and go and gi 
are independent of density. This is reminiscent of the saturation mechanism 
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Fig. 2. The equation-of-state for isospin-symmetric nuclear matter. 



in the Brueckner approach, which, to a large extent, is due to the second- 
order tensor contribution [1]. We note however that there is no one-to-one 
correspondence between the pion-exchange diagrams in the two approaches. 
For instance, the short-range part of the iterated one-pion-exchange is in our 
approach, through the renormalization procedure, subsumed in the zero-range 
interaction. Since the coupling constants go and g± are adjusted to reproduce 
the properties of nuclear matter, this part of the one-pion-exchange interac- 
tion is effectively treated non-perturbatively. In the Brueckner approach, on 
the other hand, the divergencies are regularized by form factors. Consequently, 
there the corresponding finite contributions to the energy are associated with 
the original pion-exchange diagrams. Since the one-pion-exchange is pertur- 
bative only once the short-range part has been removed, it must be iterated 
to all orders in the Brueckner approach. 

Our equation-of-state is useful for applications, where the pion dynamics plays 
an important role. As an example, we consider the quark condensate in nu- 
clear matter. The quark condensate, (qq), is an order parameter for the spon- 
taneously broken chiral symmetry. Thus, its dependence on baryon density, 
indicates to what extent chiral symmetry is restored in nuclear matter. Fur- 
thermore, it is an important ingredient in the QCD sum rules [18] and in the 
Brown- Rho scaling approach [19]. According to the Feynman-Hellman the- 
orem the quark condensate can be extracted unambiguously from the total 
energy per particle, + E(p), of nuclear matter, once its dependence on the 
current quark mass is known: 
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(qq)(p)-(qq)(0)=p- (m N + E(p, m Q )) . (17) 

OITLq 

Here mg = (m u + m d )/2 is the average of the u and d quark masses. In 
our approach the contribution of the pion degrees of freedom to the energy is 
included in a chirally consistent manner to order Q 4 . This allows us to estimate 
the effect of correlations on the quark condensate in nuclear matter to leading 
and subleading order. 

It is convenient to consider the relative change of the quark condensate 



{qq)(p) _ 1 z wNP a n (p)p 

<??>(0) mlf* 2m«tt' { ] 

The penultimate term in (18) is the modification of the quark condensate in 
a free Fermi gas of nucleons, where 



E wJV = m Q (N\qq\N) = m Q ^ (19) 
' dmQ 

is the pion-nucleon sigma term. Finally, the last term in (18) is due to the 
dependence of the interaction energy on the pion mass 



a ^- 2 <H§)dk E[p - mQ) + ° ^ <k E{p - m '» (20) 

where we use the Gell-Mann-Oakes-Renner relation to convert the depen- 
dence on the current quark mass tjiq to a dependence on the pion mass. We 
emphasize that the second term on the right-hand side of (18), which was first 
written down in refs. [20-22], does not probe the interactions in nuclear mat- 
ter. Therefore the linear density dependence of the quark condensate, which 
results when only this term is retained, should be considered with caution. It 
is a priori not clear that this term is the most important one at the saturation 
density and beyond. 

Note that here we neglect any implicit m n dependence of the effective couplings 
go and gi, since naive counting arguments suggest that the corresponding con- 
tribution to the energy per particle would be of order Q 5 . In Fig. 3 we show 
the resulting quark condensate in nuclear matter. We confront the 'leading' 
term, which is proportional to the sigma term S^jv — 45 MeV, with the full 
result given by (18). The inclusion of pionic interaction effects counteracts 
the reduction of the condensate due to the leading term. Our results confirm 
calculations performed within the Brueckner [23] and Dirac-Brueckner [24] ap- 
proach qualitatively insofar that the nuclear many-body system reacts against 
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Fig. 3. The quark condensate in isospin symmetric nuclear matter. 

chiral symmetry restoration. This result has important implications for the 
restoration of chiral symmetry in matter at finite baryon density. 
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